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Abstract: Online y-spectrometry systems for inland waters, most of which extract samples in situ and in real time, 
are able to produce reliable activity concentration measurements for waterborne radionuclides only when they are 
distributed relatively uniformly and enter into a steady-state (SS) diffusion regime in the measurement chamber. To 
protect residents’ health and ensure the safety of the living environment, better timeliness is required for this 
measurement method. To address this issue, this study established a mathematical model of the online water 
y-spectrometry system so that rapid warning and activity estimates can be obtained for water under 
non-steady-state (NSS) conditions. In addition, the detection efficiency of the detector for radionuclides during the 
NSS diffusion process was determined by applying the computational fluid dynamics (CFD) technique in 
conjunction with Monte Carlo simulations. On this basis, a method was developed that allowed the online water 
-spectrometry system to provide rapid warning and activity concentration estimates for radionuclides in water. 
Subsequent analysis of the NSS-mode measurements of “°K radioactive solutions with different activity 
concentrations determined the optimum warning threshold and measurement time for producing accurate activity 
concentration estimates for radionuclides. The experimental results show that the proposed NSS measurement 
method is able to give warning and yield accurate activity concentration estimates for radionuclides 55.42 and 
69.42 min after the entry of a 10 Bq/L *°K radioactive solution into the measurement chamber, respectively. These 
times are much shorter than the 90 min required by the conventional measurement method. Furthermore, the NSS 
measurement method allows the measurement system to give rapid (within approximately 15 min) warning when 
the activity concentrations of some radionuclides reach their respective limits stipulated in the Guidelines for 
Drinking-water Quality of the WHO, suggesting that this method considerably enhances the warning capacity of in 


situ online water y-spectrometry systems. 
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1. Introduction 

With the growing scale of nuclear energy and the extensive application of non-power nuclear 
technology, radioactive effluents discharged from nuclear facilities increase the risk of radioactive 
contamination of receiving waters. The radioactivity levels of these effluents have therefore 
become a focus of efforts to monitor the radiation environments near nuclear facilities.!!*! The 
main radionuclides measured by water y-spectrometry include natural radionuclides (e.g., the U, 
Th, and Ac series as well as *°K) and artificial radionuclides (e.g., 8Co, Co, ?!I, °4Cs, 137Cs, and 
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110m 4 g),[4+6] Two methods, laboratory analysis of samples collected at regular intervals 7!!! and in 


situ water radioactivity measurements |? 131, 


are primarily employed to monitor radioactive 
effluents discharged from nuclear facilities. In situ water radioactivity measurement is a widely 
used method owing to its advantages over laboratory analysis, such as a short measurement period, 
relatively simple procedure, and real-time monitoring. Based on the measurement environment, 


online water y-spectrometry systems can be categorized into marine and inland-water radioactivity 


measurement systems. 

Online water y-spectrometry systems for measuring marine radioactivity can be 
approximately classified into float-, cruise-, and pull-type systems. Float-type online water 
y-spectrometry systems include the MARNET marine y-radiation monitoring network developed 
by the German Federal Ministry of the Environment,!!*! the POSEIDON floating monitoring 
system developed by the Hellenic Center for Marine Research (HCMR),|!> '*! the Sentinel-type 
marine radioactivity monitoring instrument jointly developed by the Marine Environment 
Laboratory of the International Laboratory of Marine Radioactivity (IAEA-MEL) and the 
Radiological Protection Institute of Ireland,''7! and the KATERINA marine radioactivity detection 
system developed by the HCMR.!'*! Cruise-type systems include the underwater radiation spectral 
identification system developed by the United States Department of Energy’s Remote Sensing 
Laboratory!'?! and the REM-10 series underwater radiation control system developed by the 
Kurchatov Institute of the Russian Research Center.'*° *! Pull-type systems include the dual-probe 
marine in situ y-radiation monitoring system developed by the IAEA-MEL!* >l and the 
Radiometric Environment Survey and Quantification (RESQ) towed y-ray spectrometer developed 
by the University of Tokyo.!*! Monitoring networks have been established for some online water 
y-spectrometry systems to monitor radioactivity in key waters. Marine radioactivity measurement 
systems are generally unequipped with effective shielding materials owing to the restrictions 
imposed by the measurement conditions. The relatively high ambient radiation in the background 
requires a continuous measurement over more than 24 h to reach the lower detection limit of a 
marine radioactivity measurement system.|®: >: 26] 

Inland water is susceptible to significant disturbances and contains a large amount of 
sediment. Consequently, online water-spectrometry systems developed specifically for inland 
waters mostly extract samples in situ and in real time’’*?!, which are subsequently filtered and 
directed to a measurement chamber for radioactivity measurements. Waterborne radionuclides 
follow a non-steady-state (NSS) diffusion regime (i.e., a non-uniform distribution pattern) during 
the initial period after entering the measurement chamber. Because an online water-spectrometry 
system cannot accurately quantify radionuclides under NSS conditions, a wait time must be set to 
allow the radionuclides to transition from the NSS to a steady-state (SS) diffusion pattern (i.e., a 
uniform distribution pattern), after which their activity concentrations can be accurately measured. 

Conventional quantitative analysis methods quantify radionuclides based on the stable 
detection efficiencies for radionuclides in an SS diffusion pattern. Generally, the Monte Carlo 
simulation technique and validation measurements of standard radionuclide solutions are used to 
establish the detection efficiency curves for characteristic y-rays at different energy levels as well 
as to calculate the minimum detectable activity concentration, based on which radionuclides in 
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water are quantified.| 
waters that extracts samples in situ and in real time, the waiting time for the radionuclides in a 
water sample inside the measurement chamber to transition from an NSS to an SS distribution 
pattern directly increases the time required by the system to estimate the radionuclide activity 
concentrations in the water sample. This reduces the system timeliness in monitoring and giving 
warning of water radioactivity. Consequently, this study presents a method that enables online 
water-spectrometry systems designed for inland waters to provide rapid warning and activity 
estimates under NSS conditions, thus being able to protect the health of residents and ensure the 


safety of the living environment in areas near nuclear facilities. 


2. Instruments and methods 
2.1 Online water y-spectrometry system 

The online water y-spectrometry system used in this study extracts samples in situ and in real 
time. Figure 1 shows the structure of this measurement system, which consists primarily of a 
CeBr; detector, measurement chamber (made of low-background lead and oxygen-free copper), 
multichannel analysis unit, data storage system, communication unit, continuous water sampling 
device, and pollution tank. The water sample is extracted from the monitored waters and directed 
into the measurement chamber at a rate of 0.12 L/s through the continuous sampling device 
equipped with a prefiltration function. Both the inlet and outlet tubes of the measurement chamber 
have a diameter of 20 mm. Under the actual measurement conditions, the water in the 
measurement chamber is in continuous flow, the monitored waters enter through the inlet tube at 
the top and exit through the outlet tube at the bottom, and the measurement chamber is filled with 
water during the measurement. 
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Figure 1 Main structure of the online water y-spectrometry system 


The measurement unit includes the CeBr, detector and measurement chamber. The CeBr3 
detector, which is composed of three CeBr; scintillators, each with dimensions of 45 mm x 50 
mm as shown in Figure Sl(a) in the supplementary material, is used to improve the detection 
efficiency and responds to energy levels ranging from 30 keV to 3 MeV with a relative detection 
efficiency higher than 132% (@3 in Nal(TI), 1.33 MeV, Co). The measurement chamber has 
external dimensions of 740 mm x 720 mm and internal dimensions of 9496 mm x 596 mm. The 
measurement chamber has a supporting outer layer made of stainless steel and a 100-mm-thick 
internal layer made of lead. The inner side of the lead layer is covered with a 2-mm-thick 
oxygen-free copper layer and an ultrathin impermeable layer. The measurement chamber has a 
volume of 115 L. Figure S1(b) shows its detailed structure. The energy resolution of the 
measurement unit at the characteristic energy (662 keV) of '3’Cs is better than 4.5%. The lower 
detection limit for the 12 h measurement is 0.2 Bq/L for !3’Cs. 


2.2 Non-steady-state (NSS) measurement model 

When a water sample inside the measurement chamber follows the SS diffusion regime, the 
radionuclides contained in the sample are distributed uniformly inside the measurement chamber, 
with activity concentrations identical to those of the external source waters. Under these 
conditions, the detection efficiency of the measurement system for the characteristic y-rays 
emitted from a certain radionuclide in the water is constant. In the SS measurement mode, the full 


energy peak (FEP) count rate for y-rays can be described using the following equation: 


n,,(t)=C-V-e, -P i 


where nss(t) is the FEP count rate for the characteristic y-rays emitted from the radionuclide at time 
t, C is the activity concentration of the radionuclide in the measurement chamber (Bq/m°), V is the 
volume of the measurement chamber (m°), £y is the detection efficiency for the characteristic 
y-rays emitted from the radionuclide in the measurement chamber, t is the measurement time (s), 
and P, is the emission probability of the characteristic y-rays emitted from the radionuclide. 

However, the distribution of the radionuclides inside the measurement chamber of an online 
water y-spectrometry system in the NSS measurement mode changes dynamically. Consequently, 
the detection efficiency of the measurement system for the characteristic y-rays emitted from a 
certain radionuclide in the water also changes dynamically. Therefore, when a measurement 
system is in NSS measurement mode, the FEP count rate for the characteristic y-rays emitted from 
a radionuclide can be described using the following equation: 


ys (t)= A(t): e(¢)-P (2) 


where Mnss(t) is the FEP count rate for the y-rays at time t, A(t) is the activity of the radionuclide 
inside the measurement chamber at time t (Bq), and e(t) is the detection efficiency for the 
characteristic y-rays emitted from the radionuclide inside the measurement chamber at time t. 

For an online water y-spectrometry system in NSS measurement mode, because its detection 
efficiency for the characteristic y-rays emitted from a certain radionuclide in the water is not 
constant, its FEP count for the y-rays over a period of time from the starting moment to t, M(t), can 
be obtained through an integration of Equation (2): 


N(t)=J n(ede= | A(t) e(t): Pat (3) 


In the NSS measurement mode, the radionuclide in water that enters the measurement 
chamber does not completely remain inside the measurement chamber but is instead partially 
emitted from the measurement chamber with the exiting water. The calculation of the activity A(t) 
of a radionuclide in the measurement chamber requires the introduction of the proportion that 
remains in the measurement chamber, f(t), into Equation (3). f(t) represents the ratio of the activity 
of the radionuclide inside the measurement chamber to the total activity of the cumulative amount 
of radionuclide that has entered the measurement chamber from the starting moment to ¢ and is 
related to the structure of the measurement chamber and water flow rate. The activity of the 
amount of the radionuclide that remains inside the measurement chamber at time ¢ can be 
calculated based on the flow rate of the radionuclide-containing water in the inlet tube of the 
measurement chamber of the online water y-spectrometry system as well as the diameter of the 
inlet tube using the following equation: 


A(t)=C-v-t- f(t) (4) 


where d is the diameter of the inlet tube of the measurement chamber of the online 
water-spectrometry system (m), C is the activity concentration of the radionuclide entering the 
measurement chamber (Bq/m?), and v is the flow rate of water in the inlet tube of the measurement 
chamber (L/s). 

Substituting Equation (4) into Equation (3), the FEP count for the y-rays at time t, M(t), is 
obtained as follows: 


N(t)=C-v-P it f(t)-2(t)ae (5) 


Two arbitrary time points, t; and t2, during the NSS measurement process were selected for 
the analysis and substituted into Equation (5). The activity concentration of the radionuclide in the 
water flowing into the measurement chamber was then calculated using the following equation: 


N(t)-N() 
P, v- f telt) f (at 


(6) 
The use of Equation (6) to calculate the activity concentration of a radionuclide during the 


NSS measurement process requires obtaining the variables f(t) and e(t) first. 


3. Results 
3.1 Residual percentage f(t) of radionuclide 

An analysis based on the diameter and water flow rate in the inlet tube of the measurement 
chamber revealed that the water that enters this chamber follows a turbulent flow regime. This 
presents challenges in determining the distribution patterns of the radionuclide at different time 
points and the variation in its concentration with time through mathematical description and direct 
measurements. Therefore, computational fluid dynamics (CFD)°” was used to simulate the 
distribution pattern of radionuclides in water after entry into the measurement chamber. Based on 
this, f(t) in the NSS measurement mode was calculated. Figure 2 shows a schematic of the 
geometry of the input model constructed in the CFD simulation based on the actual dimensions of 
the measurement chamber of the online water y-spectrometry system. The measurement chamber 
was composed of a cylinder with dimensions of 9496 mm x 596 mm and a hollow cylinder with 
dimensions of 140 mm x 363 mm. The hollow cylinder acted as a cavity to house the CeBr3 
detector. Both the inlet and outlet tubes of the measurement chamber had a diameter of 20 mm. 


The water flow rates in the two tubes were 0.12 L/s. 
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Figure 2 Schematic of the geometry of the input model of the measurement chamber 


Figure 3(a) shows the analysis of the average radionuclide concentration across the 
measurement chamber during the diffusion process for the 100 Bq/L “°K radioactive solution 
based on the domain probe module in the CFD simulation, and the average radionuclide 
concentration does not consider the effects of radionuclide decay. The radionuclide concentration 
in the water inside the measurement chamber becomes consistent with that in the source water 


only after a diffusion time of 90 min, suggesting that a wait time of at least 90 min is required to 


accurately quantify the radionuclide activity concentration in the water using the online water 
y-spectrometry system coupled with the conventional measurement method. 
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(a) and (b). 
Figure 3 Simulated diffusion of the radionuclides in the measurement chamber (curves depicting the variation in 
the (a) average concentration and (b) f(t) of the radionuclides inside the measurement chamber with the diffusion 
time) 


Based on the analysis of the variation in the average radionuclide concentration inside the 
measurement chamber with the diffusion time, f(t) can be calculated using the following equation 
for the online water y-spectrometry system in the NSS measurement mode: 

fy- OF 
vet-C (7) 
where v (= 0.12 L/s) is the water flow rate in the inlet and outlet tubes during the dynamic 
measurement process, C (= 100 Bq/L) is the radionuclide concentration in the source water outside 
the measurement chamber, C(t) is the average radionuclide concentration inside the measurement 
chamber at time ¢ (Bq/L), and V (= 0.115 m’) is the volume of the measurement chamber. 

Figure 3(b) shows the ft) values at different diffusion times, calculated using Equation (7). 
An analysis of the calculation results provided a fitting equation for f(t) and the diffusion time t 
(R? = 0.998), as shown in Equation (8). 

-t 
f (t) = 0.1343 + 0.6594 x 2077910" (8) 
3.2 Monte Carlo simulations of the detection efficiency 

To calculate the time-varying detection efficiency of the detector inside the measurement 
chamber for the characteristic y-rays emitted from a radionuclide in the water in NSS 
measurement mode, a physical model of the source distribution corresponding to the radionuclide 
inside the measurement chamber was constructed based on the data obtained in Section 3.1 for the 
concentration distribution of the radionuclide during the NSS diffusion process. The concentration 
distribution of the radionuclides inside the measurement chamber during the NSS diffusion 
process cannot be expressed as a function. Thus, the measurement chamber was divided into a 
finite number of water elements, which were then filled with the radionuclide concentrations in the 
water at different time points to simulate the detection efficiency of the measurement system under 
NSS conditions using the Monte Carlo technique in the open-source software Geant4. To ensure 
the accuracy of the data for the distribution of the radionuclide concentration inside the 
measurement chamber and considering the computational efficiency, the measurement chamber 


was divided into 17 layers, as shown in Figure S2(a). Subsequently, as shown in Figure S2(b), the 


water in each layer was divided into multiple cuboids/fan-shaped curved bodies. The entire 
measurement chamber was divided into 1864 water elements. 

The radionuclide concentration of each water element in the measurement chamber was 
calculated using the water elements employed as domain probes in the CFD simulation. Figure 4 
shows the distribution patterns of the radionuclide concentrations in each water element in the 7th 
layer of the measurement chamber 5, 10, 30, and 60 min after the radionuclide-containing water 
entered the measurement chamber. The radionuclide concentration in each water element 
continuously increased with radionuclide diffusion time. Similar radionuclide concentrations were 
reached in the different water elements in this layer 60 min after the radionuclide entered the 
measurement chamber. The mass of the radionuclides at different locations inside the 
measurement chamber was determined by multiplying the radionuclide concentrations obtained 
for all water elements by the volume of the corresponding area. The distribution patterns of the 
radionuclide concentration inside the measurement chamber at different diffusion times were then 


determined by normalizing the results obtained for all the water elements. 


(O) (d) 
Figure 4 Distribution patterns of the radionuclide concentration in each domain probe area in the 7th layer of the 
measurement chamber at diffusion times of (a) 5, (b) 10, (c) 30, and (d) 60 min 


The geometry of the physical model constructed to calculate the detection efficiency through 
Monte Carlo simulation agreed completely with that of the measurement chamber. The 
radionuclide-containing water inside the measurement chamber was divided into the same number 
of water elements as those shown in Figure S2 and sampled as the source term. The radionuclide 
concentration in each water element was set based on the corresponding normalized concentration 
determined by the CFD simulation. Figure S3 shows the physical structure of the model used to 
calculate the detection efficiency inside the measurement chamber through the simulation and 


partitioning of the model into water elements. 


Monte Carlo numerical simulations yield curves depicting the variation in the detection 
efficiencies of the system for the characteristic y-rays emitted from “°K (1460.9 keV), >8Co 
(810.77 keV), ©Co (1332.51 keV), P! (364.49 keV), '°4Cs(604.72 keV), and '°’Cs(661.66 keV) 
inside the measurement chamber with the diffusion time, as shown in Figure 5. The detection 
efficiency of the measurement system for radionuclides increased rapidly as the radionuclides 
continuously entered the measurement chamber during the initial stage. When the diffusion time 
increased beyond 900 s, the distribution pattern of each radionuclide inside the measurement 
chamber remained almost unchanged, whereas the concentration of each radionuclide 
continuously increased. Under this condition, the detection efficiency of the measurement system 


for the characteristic y-rays emitted from each radionuclide basically tends to stabilize. 
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Figure 5 Curves of the variation in the detection efficiency of the measurement system for the characteristic y-rays 
emitted from (a) K (1460.9 keV), (b) 5Co (810.77 keV), (c) Co (1332.51 keV), (d) B'I (364.49 keV), (e) “Cs 
(604.72 keV), and (f) "Cs (661.66 keV) with the diffusion time 


The detection efficiencies of the measurement system for y-rays at different energy levels in 
the NSS measurement mode obtained from Monte Carlo simulations reveal that they tend to level 
off 900 s after the entry of radionuclides into the measurement chamber. When the activity 
concentration of a radionuclide in the water is calculated, if the values of t; and tz in Equation (6) 
are both set to be greater than 900 s, the detection efficiency e(t) can be considered as a constant, 
as presented in Table 1. Under these conditions, the cumulative count and activity concentration of 
the radionuclides in the NSS measurement mode can be calculated using the following equations: 

t 
N(t)=C-v-P-e,-[ t- f(t) dt O) 
N (t) -N (t ) 


P, v-e, f t-f (e)dt 155 


C= 


Table 1 Detection efficiencies for y-rays at different energy levels in non-steady-state (NSS) measurement mode 
y-ray Energy Emission Probablity 
(keV) (%) 


Radionuclide Detection Efficiency (X 10°) 


K 1460.82 10.67 1.05 
58Co 810.71 99.45 1.39 
Co 1332.51 99.98 1.10 
maij 364.49 81.24 1.93 
134Cs 604.72 97.60 1.59 
BICs 661.66 84.62 1.53 


3.3 Experimental validation 

The relative deviation between the cumulative count of the radionuclides given by the NSS 
measurement model and the cumulative count measurement directly reflects the measurement 
accuracy of the model for the radionuclide concentration. To prevent radioactive contamination of 
the surrounding environment, *°K, an easily accessible natural radionuclide, was used in the 
validation experiment. Different masses of KCI powder were dissolved in water to prepare “°K 
radioactive solutions with six activity concentrations (5, 10, 20, 30, 40, and 60 Bg/L). 
Subsequently, the system was used to obtain the y-spectrometry measurements for the *°K 
radioactive solutions with different activity concentrations at different diffusion times in the NSS 
measurement mode, as shown in Figure S4. Figure 6 shows the measured net FEP counts for K 
radioactive solutions with different activity concentrations in the NSS measurement mode. 
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Figure 6 Measured net full energy peak (FEP) counts for the *°K radioactive solutions with activity concentrations 
of (a) 5, (b) 10, (c) 20, (d) 30, (e) 40, and (£) 60 Bq/L 


Analysis of the relative errors of the measured cumulative counts for the *°K radioactive 
solutions with different activity concentrations from the theoretical values revealed a basic 
agreement between the measurements and theoretical calculations of the radionuclide in water. 
The high statistical fluctuation in the measured counts for each “°K-containing radioactive solution 
in the first 20 min after its entry into the measurement chamber, owing to the relatively low 
activity concentration of K in the chamber, led to a relatively significant relative error of the 
measurement from the theoretical calculation. A continuous increase in the activity concentration 
of *°K in the measurement chamber decreased the relative error, which fell within +5% after the 
first 20 min for each solution. 


4. Discussion 
4.1 Predictive analysis of radionuclides in water in the ideal measurement condition 

To determine an optimum warning threshold and an optimum prediction time for the online 
water y-spectrometry system in the NSS measurement mode, the NSS measurements of the 10 
Baq/L “°K radioactive solution were selected for analysis. An ideal measurement condition refers to 
a situation in which the specific entry time of the radionuclide-containing water into the 
measurement chamber is known. Under normal circumstances, radioactivity measurements at a 
confidence level of 0.95 are considered reliable. Based on eight y-spectrometry measurements 
performed on the 10 Bq/L *°K radioactive solution in the NSS diffusion mode, confidence levels 
and intervals were calculated for the net FEP counts for “°K at 1460.9 keV at different diffusion 
times (as presented in Table S1 in the supplementary material). An average net FEP count of 143 + 
6.44 was determined for *°K at 1460.9 keV 35 min after the beginning of the diffusion of the °K 
radioactive solution inside the measurement chamber, with a confidence level of greater than 0.95. 
The FEP background count at 1460.9 keV at this time was 214 + 11.04. The ratio of the 
cumulative FEP count to the natural FEP background count was 1.67. The use of this ratio as a 
warning threshold for the water radioactivity measurement system is relatively reasonable. 

When the NSS measurement model in Equation (10) is used to estimate the radionuclide 
activity concentration in water, the time at which the warning threshold of the measurement 
system is reached can be selected as t;. Then, Equation (10) can be used to calculate the 
radionuclide activity concentrations at different time intervals At between t; and tz (as listed in 
Table 2), and the uncertainty of the activity concentration estimate is the synthetic uncertainty, 
which is mainly contributed by the statistical fluctuation of the cumulative FEP count. The 
calculation results in Table 2 reveal that at At > 15 min, the relative deviation of the radionuclide 
activity concentration in the water yielded by the NSS measurement model from the true value is 
less than 5%, suggesting that this model is effective at giving relative accurate estimates for 
radionuclide activity concentrations in water. Therefore, producing an accurate estimate of the 
radionuclide activity concentration in water requires measurements for an additional period of 15 


min after the warning threshold of the measurement system is reached. 


Table 2 Activity concentration estimates for “°K in the 10 Bq/L *°K radioactive solution at different time intervals 
in the ideal measurement condition. 


Activity concentration for °K 


Time interval A¢(min) Relative error (%) 
(Bq/L) 
1 5.71 £1.94 -42.90 
5 7.65 + 1.61 -23.50 
10 9.33 + 1.21 -6.70 
15 9.50 + 0.79 -5.00 
20 9.83 + 0.74 -1.70 
30 9.86 + 0.59 -1.40 


The time at which the warning threshold of the measurement system was reached varied with 
the activity concentration of ®K in the *°K radioactive solution. In addition, the use of the NSS 
measurement model to calculate the radionuclide activity concentration in water requires the 
determination of the diffusion time of the radionuclide when the warning threshold of the 


measurement system is reached, ¢7, and an interval of 15 min between ¢; and t2. Based on Equation 


(9), the following equation is obtained: 


t,=t,+900 (a1) 
N(t,)_ Joe fae 
N(t, ae ' 

(£) f Es) dt 5) 


Analysis of the measurements determined the net FEP counts for “°K at 1460.9 keV at t; and 
t2 (N(t2) and N(ti), respectively), as well as their ratio k. Substituting Equations (8) and (11) into 
Equation (12) gives the following higher-order equation with respect to ty: 


t t 


of? : t eC -eG4]= 
Dt tE- t Fte? tHe tst (13) 


Table 3 summarizes the values of coefficients D, E, F, G, H, and I in Equation (13). 


Table 3 Values of the coefficients in Equation (13) 


Coefficients Values 
D 6.7135 X 102X(1-k) 
1.2084 x 10? 
8.8851 X 10?-1.3701 X10? Xk 
2.0779 X 103 
7.9970 X 105 
5.4379 X 104 
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In summary, the following procedure can be used to calculate the radionuclide activity 
concentration in radioactively contaminated water. First, the ratio of the net FEP count for the 
characteristic y-rays emitted from the radionuclide at the time when the preset warning threshold 
of the measurement system is reached to that at 15 min after reaching the warning threshold, k, is 
calculated. Subsequently, substituting k into Equation (13) determines the diffusion time of the 
radionuclide inside the measurement chamber when the warning threshold of the measurement 
system is reached, t;. Finally, substituting t; and tz into the NSS measurement model in Equation 
(10) provides an activity concentration estimate for the radionuclide in the monitored water. 


4.2 Predictive analysis of radionuclides in water in actual measurement conditions 

In actual measurement conditions, to eliminate the impact of statistical fluctuations on 
measurements, the ratio of the cumulative FEP count for the characteristic y-rays emitted from the 
radionuclide of interest in the water to the corresponding natural background count within a fixed 
time interval T is used to determine whether the warning threshold of the online water 
y-spectrometry system for the radionuclide activity concentration is reached. In this study, a 10 
Bq/L *°K radioactive solution was selected for the measurement. Figure 7 shows the variation in 
the FEP count rate given by the online water y-spectrometry system in NSS measurement mode 
for ®K at 1460.9 keV with the measurement time. When the NSS measurement model is used to 
estimate the radionuclide activity concentration in water, accurately determining whether the 
measurement reaches the warning threshold requires the selection of a suitable time interval T to 


facilitate the acquisition of reliable measurement data. 
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Figure 7 Variation in the FEP count for *°K at 1460.9 keV with the measurement time 


As shown in Figure 7, the measurement system must perform continuous measurements of 
the monitored water under actual measurement conditions. The use of the NSS measurement 
model to estimate the radionuclide activity concentration in water assumes that the radionuclide 
enters the measurement chamber at time ¢ = 0 and that the warning threshold of the measurement 
system is reached at time ¢t;. The following can be obtained based on the warning threshold 
selected for ideal measurement conditions. 


NUJEN =1.67 


j (14) 
where M(t) is the net FEP count for ®K at 1460.9 keV after the radionuclide diffuses inside the 
measurement chamber for a period from the starting moment (t = 0) to t; and Ng is the FEP 
background count for “°K at 1460.9 keV within the measurement time T. 

To select a suitable measurement time 7, Equation (15) is used to calculate the measurement 
accuracies for the FEP background counts at 1460.9 keV on the background y-spectrometry given 


by the measurement system within different measurement times. 


%)= a x 100% 

n(%) Jn T 
where 7 is the measurement accuracy and N is the FEP background counts at 1460.9 keV within 
the given time. High-accuracy measurements require a value of 7 less than 5%. 

Table S2 lists the calculated 7 values for the FEP background counts at 1460.9 keV on the 
background y-spectrometry given by the measurement system within different measurement times. 
At T= 65 min, an 7 value below 5%, as required for high-accuracy measurements, can be reached. 
Analysis of the 11 background measurements taken by the measurement system reveals that the 
mean Nz at 1460.9 keV at T = 65 min is 417.06, with a mean squared error o of 8.98. Therefore, 
when M(t) > 30, the confidence level can reach above 99.7% when ft; is set as the warning time. 

The net FEP count for K at 1460.9 keV, N(tı), and the diffusion time of the radionuclide 
inside the measurement chamber, tı, at T = 65 min can be calculated for the 10 Bq/L K 
radioactive solution using Equations (14) and (9), respectively: 


et = 279.43 


t, = 3325 ie 


Under this condition, the calculated value of N(tı) is much higher than 30. The preset 
warning threshold allows for effective warning of radioactive contamination in water. In addition, 
for the 10 Bq/L *°K radioactive solution, the measurement system gives warning 55.42 min after 
the entry of the radionuclide into the measurement chamber. Subsequently, Equation (10) is used 
to calculate the radionuclide activity concentrations in the water at different time intervals between 
tı and t2, At (as listed in Table S3). When the measurement was continued for an additional period 
of 14 min after the warning threshold of the measurement system was reached, the relative 
deviation of the radionuclide activity concentration provided by the NSS measurement model 
from the true value was less than 5%. 

Based on the above findings, the measurement time T and At should be set to 65 and 14 min, 
respectively, in the NSS measurement model established in this study for the selected online water 
y-spectrometry system when used in actual measurement conditions. The following procedure was 
used to calculate the radionuclide activity concentration in the monitored water: The time at which 
the warning threshold of the measurement system was reached was defined as t7. In addition, the 
net FEP count for the characteristic y-rays at t;, N(ti), and that at the time ¢2 after 14 min of 
continual measurements, N(t2), were recorded. Then, the N(t2)/N(tı) ratio, k, was substituted into 
Equation (13) to calculate the diffusion time of the radionuclide inside the measurement chamber 
i.e., t1) when the warning threshold of the measurement system was reached at the At of 14 min. 
Finally, substituting the obtained calculation results into Equation (10) provided an activity 


concentration estimate for the radionuclides in the monitored water. 


Table 4 Warning times and relative deviations of the estimates given by the NSS measurement model from the true 
values for “°K radioactive solutions with different activity concentrations 


True activity concentration Warning times Estimated activity concentration Relative error 
(Bq/L) (min) (Bq/L) (%) 
10 55.42 9.61 + 0.71 -3.90 
20 34.39 20.31 1.13 1.53 
30 27.54 29.46 + 1.46 -1.79 
40 20.97 39.26 + 1.85 -1.85 
60 18.33 59.49 + 2.56 -0.85 
80 15.02 78.63 + 3.85 -1.71 
100 14.00 99.11 + 3.87 -0.89 


Table 4 summarizes the times required by the measurement system to reach the warning 
threshold at actual measurement conditions (T = 65 min; At= 14 min) for “°K radioactive solutions 
with different activity concentrations as well as the relative deviations of the estimated activity 
concentrations of the radionuclide in the water from the true values. The measurements show that 
the NSS measurement model established in this study is effective in accurately estimating the 
radionuclide activity concentration in water under actual measurement conditions. 

Further analysis of the detection efficiencies of the measurement system for characteristic 
y-rays at different energy levels simulated using the Monte Carlo method reveals that their ratios 
are fixed when they level off. Therefore, the activity concentrations of different radionuclides at 
the same count rate can be mutually converted based on the emission probability and detection 
efficiency of the measurement system for their characteristic y-rays at different energy levels. The 


analysis based on the fitting equation for the variation in the warning time with an activity 


concentration of K in water in Figure S5 shows that the measurement system requires 
approximately 15 min to provide a warning when the activity concentration of each radionuclide 
in Table 5 reaches the limit stipulated in the Guidelines for Drinking-water Quality (GDWQ) of 
the World Health Organization.!**! 


Table 5 Times required by the measurement system to give warning when the activity concentrations of different 
radionuclides reach their respective limits stipulated in the Guidelines for Drinking-water Quality (GDWQ) of the 
World Health Organization. 


Corresponding “°K 


y-ray Energy Guidance level Early warning time 


Radionuclide (keV) (Ba/L) a ni (iain 
Co 810.71 100 1233.86 15.00 + 0.13 
Co 1332.51 100 981.64 15.00 + 0.11 

BIT 364.49 10 139.95 15.01 +0.17 
134Cs 604.72 10 138.51 15.01 + 0.12 
Cs 661.66 10 115.56 15.05 + 0.14 


5. Conclusions 

This study presents an NSS measurement model for radionuclide activity concentrations in 
water to complement online water y-spectrometry systems designed for inland waters that mostly 
extract samples in situ and in real time. On this basis, a method is developed to facilitate the online 
water y-spectrometry system to give rapid warning and activity estimates in NSS conditions. The 
experimental results show that when coupled with the NSS measurement model, the measurement 
system is able to give warning and obtain an accurate estimate for the activity concentration of the 
radionuclide (K) in water 55.42 and 69.42 min, respectively, after the entry of a 10 Bq/L “K 
radioactive solution (corresponding to a 0.87 Bq/L '3’Cs radioactive solution) into the 
measurement chamber, and the relative error of activity concentration is below + 5%. These times 
are much shorter than the 90 min required for the conventional measurement method. In addition, 
the NSS method allows the measurement system to provide rapid (within approximately 15 min) 
warnings when the activity concentrations of some radionuclides reach their respective limits, as 
stipulated in the GDWQ of the World Health Organization. These results show that the method 
developed in this study allows online water y-spectrometry systems to give rapid warning and 
produce activity estimates in NSS conditions. Therefore, this method significantly improves the 
timeliness of the monitoring of the radioactivity levels of inland waters that receive effluent from 
nearby nuclear facilities. It can provide rapid and reliable measurement data to facilitate daily 


monitoring of the radioactivity levels of inland waters and monitoring of nuclear emergencies. 
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